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Abstract. - Magnetic fields and magnetic impurities are each known to suppress super- 
conductivity. However, as the field quenches (i.e. polarizes) the impurities, rich consequences, 
including field-enhanced superconductivity, can emerge when both effects are present. For the 
case of superconducting wires and thin films, this field-spin interplay is investigated via the 
Eilenberger-Usadel scheme. Non-monotonic dependence of the critical current on the field (and 
therefore field-enhanced superconductivity) is found to be possible, even in parameter regimes 
in which the critical temperature decreases monotonically with increasing field. The present 
work complements that of Kharitonov and Feigel'man, which predicts non-monotonic behavior 
of the critical temperature. 



Introduction. - In their classic work, Abrikosov and Gor'kov [1] predicted that unpolar- 
ized, uncorrelated magnetic impurities suppress of superconductivity, due to the de-pairing ef- 
fects associated with the spin-exchange scattering of electrons by magnetic impurities. Among 
their results is the reduction, with increasing magnetic impurity concentration, of the super- 
conducting critical temperature T c , along with the possibility of "gapless" superconductivity 
in an intermediate regime of magnetic-impurity concentrations. The latter regime is realized 
when the concentration of the impurities is large enough to eliminate the gap but not large 
enough to destroy superconductivity altogether. Not long after the work of Abrikosov and 
Gor'kov, it was recognized that other de-pairing mechanisms, such as those involving the cou- 
pling of the orbital and spin degrees of freedom of the electrons to a magnetic field, can lead 
to equivalent suppressions of superconductivity, including gapless regimes [2-5] . 

Conventional wisdom holds that magnetic fields and magnetic moments each tend to sup- 
press superconductivity (see, e.g., Ref. [6]). Therefore, it seems natural to suspect that any 
increase in a magnetic field, applied to a superconductor containing magnetic impurities, would 
lead to additional suppression of the superconductivity. However, very recently, Kharitonov 
and Feigel'man [7] have predicted the existence of a regime in which, by contrast, an increase 
in the magnetic field applied to a superconductor containing magnetic impurities leads to 
a critical temperature that first increases with magnetic field, but eventually behaves more 
conventionally, decreasing with the magnetic field and ultimately vanishing at a critical value 
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of the field. Even more strikingly, they have predicted that, over a certain range of concen- 
trations of magnetic impurities, a magnetic field can actually induce superconductivity out of 
the normal state. 

The Kharitonov-Feigel'man treatment focuses on determining the critical temperature by 
determining the linear instability of the normal state. The purpose of the present Letter is 
to address properties of the superconducting state itself, most notably the critical current 
and its dependence on temperature and the externally applied magnetic field. The approach 
that we shall take is to derive the (transport-like) Eilenberger-Usadel equations [9,10], by 
starting from the Gor'kov equations. We account for the following effects: potential and spin- 
orbit scattering of electrons from non-magnetic impurities, and spin-exchange scattering from 
magnetic impurities, along with orbital and Zceman effects of the magnetic field. In addition 
to obtaining the critical current, we shall recover the Kharitonov-Feigel'man prediction for the 
critical temperature, as well as the dependence of the order parameter on temperature and 
applied field. In particular, we shall show that not only are there reasonable parameter regimes 
in which both the critical current and the transition temperature vary non-monotonically 
with increasing magnetic field, but also there are reasonable parameter regimes in which only 
the low-temperature critical current is non-monotonic even though the critical temperature 
behaves monotonically with field. The present theory can be used to explain certain recent 
experiments on superconducting wires [8]. 

Before describing the technical development, we pause to give a physical picture of the 
relevant de-pairing mechanisms. First, consider the effects of magnetic impurities. These 
cause spin-exchange scattering of the electrons (including both spin-flip and non-spin-flip 
terms, relative to a given spin quantization axis), and therefore lead to the breaking of Cooper 
pairs [1]. Now consider the effects of magnetic fields. The vector potential (due to the 
applied field) scrambles the relative phases of the partners of a Cooper pair, as they move 
diffusively in the presence of impurity scattering (viz. the orbital effect), which suppresses 
superconductivity [2,3]. On the other hand, the field polarizes the magnetic impurity spins, 
which decreases the rate of exchange scattering (because the spin-flip term is suppressed more 
strongly than the non-spin- flip term is enhanced), thus diminishing this contribution to de- 
pairing [7] . In addition, the Zecman effect associated with the effective field (coming from the 
applied field and the impurity spins) splits the energy of the up and down spins in the Cooper 
pair, thus tending to suppress superconductivity [6] . We note that strong spin-orbit scattering 
tends to weaken the de-pairing caused by the Zeeman effect [5] . Thus we see that the magnetic 
field produces competing tendencies: it causes de-pairing via the orbital and Zeeman effects, 
but it mollifies the de-pairing caused by magnetic impurities. This competition can manifest 
itself through the non-monotonic behavior of observables such as the critical temperature and 
critical current. In order for the manifestation to be observable, the magnetic field needs to be 
present throughout the samples, the scenario being readily accessible in wires and thin films. 

The model. - We take for the impurity-free part of the Hamiltonian the BCS mean-field 
form [6,11]: 

H = ~f rfr ^«( V -7 A ) 2 ^ + Y / dr ((^J>^a+^J<^a>)-/i J dr^ a , 

(1) 

where ^>£(r) creates an electron having mass m, charge e, position r and spin projection a, 
A is the vector potential, c is the speed of light, /i is the chemical potential, and Vq is the 
pairing interaction. Throughout this Letter we shall put % = 1 and ks = 1- Assuming the 
superconducting pairing is spin-singlet, we may introduce the complex order parameter A, 
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via 

-WaV/3> = KpA, VoiMp) = ia y a/3 A*, (2) 

where <r^' z are the Pauli matrices. We assume that the electrons undergo potential and 
spin-exchange scattering from the magnetic impurities located at a set of random positions 
{xi}, in addition to undergoing spin-orbit scattering from an independent set of impurities 
or defects located at an independent set of random positions {yj}, as well as being Zeeman 
coupled to the applied field: 

Hint = J driliV^p, (3a) 

with V a p being given by 

V a p = {u 1 (r-x i )5 a p + u 2 (r-x l )Si-a a p}+^^{Vv so (r-y j )-(a a pxp)}+[i B B<Jap, (3b) 

where Si is the spin of the z-th magnetic impurity and where, for simplicity, we have attributed 
the potential scattering solely to the magnetic impurities. We could have included potential 
scattering from the spin-orbit scattering centers, as well as potential scattering from a third, 
independent set of impurities. However, to do so would not change our conclusions, save for 
the simple rescaling of the mean-free time. We note that cross terms, i.e. those involving 
distinct interactions, can be ignored when evaluating self-energy [5,7]. Furthermore, we shall 
assume that the Kondo temperature is much lower than the temperature we are interested in. 

The impurity spins interact with the applied magnetic field through their own Zeeman 
term: 

H z = -u s S z (4) 

where oj s = gsHsB, and g s is the impurity-spin g-factor. Thus, the impurity spins are not 
treated as static but rather have their own dynamics, induced by the applied magnetic field. 
We shall approximate the dynamics of the impurity spins as being governed solely by the 
applied field, ignoring any influence on them of the electrons. Then, as the impurity spins are 
in thermal equilibrium, we may take the Matsubara correlators for a single spin to be 

(T T S+{n)S-fo)) = TV D+re-™'^-^, (5a) 

(T T S-(n)S+(T 2 )) = TV D-+e-^'^-^, (5b) 

* ■* ii)' 



(TrS'WS'fo)) = d z = (S*) 2 , (5c) 

where w' (= 2irnT) is a bosonic Matsubara frequency, ~ denotes a thermal average, and 

D+r = + u a ), D~+ = 2S^/(+ilu' + u s ). (6) 

We shall ignore correlations between distinct impurity spins, as their effects are of the second 
order in the impurity concentration. 

To facilitate the forthcoming analysis, we define the Nambu-Gor'kov four-component 
spinor (sec, e.g., Refs. [5,12]) via 

= (^}(r,r),V{(r,T),VT(r,T),Vi(r,T)). (7) 
Then, the electron-sector Green functions are defined in the standard way via 

g,(i : 2) = - ( r T * <( i)*t (2)> s |a - 2)) , (8) 
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where G, &, F, and .F^ are each two-by-two matrices (as indicated by the * symbol), being 
the normal and anomalous Green functions, respectively. As the pairing is assumed to be 
singlet, F is off-diagonal whereas G is diagonal. 

Eilenberger-Usadel equations. The critical temperature and critical current are two 
of the most readily observable quantities. As they can be readily obtained from the Eilen- 
berger and Usadel equations, we shall focus on these formalisms. A detailed derivation will 
be presented elsewhere. The procedure is first to derive Eilenberger equations [9], and then, 
assuming the dirty limit, to obtain the Usadel equations. The self-consistency equation be- 
tween the anomalous Green function and the order parameter naturally leads, in the small 
order-parameter limit, to an equation determining the critical temperature. Moreover, solv- 
ing the resulting transport-like equations, together with the self-consistency equation, gives 
the transport current, and this, when maximized over supcrfluid velocity, yields the critical 
current. 

To implement this procedure, one first derives the equations of motion for Q (viz. the 
Gor'kov equations). By suitably subtracting these equations from one another one arrives at 
a form amenable to a semiclassical analysis, for which the rapidly and slowly varying parts 
in the Green function (corresponding to the dependence on the relative and center-of-mass 
coordinates of a Cooper pair, respectively) can be separated. Next, one treats the interaction 
Hamiltonian as an insertion in the self-energy, which leads to a new set of semi-classical 
Gor'kov equations. These equations are still too complicated to use effectively, but they can 
be simplified to the so-called Eilenberger equations [9,13-15] (at the expense of losing detailed 
information about excitations) by introducing the energy-integrated Green functions, 

g(u, k,R)= l - J di k G(u, k, R), f(u, k,R) = ^ J d£ k F(u, k, R), (9) 

and similarly for g^(iv,k,R) and p(u,k,R). Here, to is the fermionic frequency Fourier con- 
jugate to the relative time, k is the relative momentum conjugate to the relative coordinate, 
and R is the center-of-mass coordinate. (We shall consider stationary processes, so we have 
dropped any dependence on the center-of-mass time.) However, the resulting equations do not 
determine g's and f's uniquely, and they need to be supplemented by additional normalization 
conditions [9,13-15], 

g 2 + fP = g n + f\f = i, (10) 

as well as the self-consistency equation, 

A = |5l£ /12M. (11) 

In the dirty limit (i.e. wr tr -C G and Ar tr -C F), where r tr is the transport relaxation time 
(which we do not distinguish from the elastic mean- free time), the Eilenberger equations can 
be simplified further, because, in this limit, the energy-integrated Green functions are almost 
isotropic in k. This allows one to retain only the two lowest spherical harmonics (I = 0, 1), 
and to regard the I = 1 term as a small correction (i.e. \k ■ F\ -C |F|) so that we may write 

g{u,k,R) = G{u,R) + k-G(uj,R), f{u,k,R) = F(u,R) + k-F(u,R), (12) 

where k is the unit vector along k. In this nearly-isotropic setting, the normalization conditions 
simplify to 

0^ = 1-^12^21, G\ 2 = \-F2\F\^ (13) 

and the Eilenberger equations reduce to the celebrated Usadel equations [10] for Fi 2 (w,R), 
F 21 (u,R), Fl 2 (cu,R), and F^(w,R). 
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Application to thin wires and films. - Let us consider a wire (or film) not much thicker 
than the effective coherence length. In this regime, we may assume that the order parameter 
has the form A(R) — Ae mRx , where R x is the coordinate measured along the direction of the 
current (e.g. for a wire this is along its length) and u is a parameter encoding the velocity 
of the superfiow h~u/2m. Similarly, we may assume that the semiclassical anomalous Green 
functions have a similar form: 

F 12 {u,R) = F 12 (u)e iuR * , F 21 (co,R) = F 21 (u)e iuR * , (14a) 

F\ 2 (to, R) = Fl 2 ( w ) e -*»«- , F 2 \ (cj, R) = F 2 \ Me"** . (14b) 

Together with the symmetry amongst F's (i.e. F*p — — i^L and F a p = —Fp a ) we can reduce 
the four Usadel equations for Fi2, F21, -F 12 , and F 21 to one single equation: 

to + i5 B + ^— (d-+G 2 2{u-uj')) + f— + 5gii(w) + ^-G 22 (w) 

2t b *—f \ / \T B 2 / 6T so 

in which Sg = \ibR + niU 2 (0)S z , D = D(((u — 2eA/c) 2 }} with the London gauge chosen 
and ((■ • • )) defining a spatial average over the sample thickness, D = VpTt r /3 is the diffusion 
constant. The spin-exchange and spin-orbit scattering times, tb and r so , are defined via the 
Fermi surface averages 

^- = Nomn J ^K| 2 , = N Q n so Tr J ^-\v so \ 2 p 2 F \k x fcf. (16) 

Here, N is the (single-spin) density of electronic states at the Fermi surface, rij is the concen- 
tration of magnetic impurities, n so is the concentration of spin-orbit scatterers, andpp = tuvf 
is the Fermi momentum. The normalization condition then becomes 

GnM = sgnH[l - F? 2 (")} 1/2 , G 22 (uj) = sgn(w)[l - F* 2 ^)] 1 ' 2 = G* xl {u). (17) 

Furthermore, the self-consistency condition (|11J) becomes 

Aln(T C0 /T) =7rr^((A/|w|) -F 12 (u)), (18) 

in which we have exchanged the coupling constant g for Tco, i.e., the critical temperature of 
the superconductor in the absence of magnetic impurities and fields. 

In the limit of strong spin-orbit scattering (i.e. r so <C I /to and tb), the imaginary part of 
Eq. lfl5|) is simplified to 

\S B + — ImV D-+ G(u!-ui')]ReC + —lm(GC) = 0, (19a) 
2r B ^— 'W - 3r so 

and the real part is rewritten as 

wReG + JjLRe£ u , G{lu-lu') C{lu) + G*{u) D~+(uj') C*(u - w')] 

-[S B + ^Im£ w , DZ+ G(u-u')]lmC+ (f + f ) Rc(GG) = Re G, (19b) 



F 12 {lo) 
A 
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where C = F12/A, G = Gu, and the argument w is implied for all Green function, except 
where stated otherwise. Next, we take the advantage of the simplification that follows by 
restricting our attention to the weak-coupling limit, in which F\2{lo) <C 1. Then, eliminating 
G in Eq. Ijl9|) using Eqs. I|17(l . and expanding to third order in powers of F, one arrives at 
an equation for F that is readily amenable to numerical treatment. The quantitative results 
that we now draw are based on this strategy. ( 1 ) 

Results for the critical temperature. - These can be obtained in the standard way, i.e., 
by (i) setting u = and expanding Eqs. i|19|) to linear order in F (at fixed A), and (ii) setting 
A — > and applying the self-consistency condition. Step (i) yields 



M + f " + T(((^) >> + ^Hl ReC H-i-4E^r^ Re ^--')' ( 20a ) 

to' s 

where 



sM^s B -^- £ W (20b) 

Tr> *■ * 1.1' -i-f.l^ 



W e >|W |>M 



in which a cutoff uj c has been imposed on to', and 



t, s «- + L E (21) 



tb Tb . frr . u> + w 

w < a? 



This is essentially the Cooperon equation in the strong spin-orbit scattering limit, first derived 
by Kharitonov and Feigel'man [7], up to an inconsequential renormalization of 8b- 
Step (ii) involves solving the implicit equation 



^ = ^£[^4(c ( .) + c>) 



(22) 



the solution of which is T = Tq- 

Figure ^ shows the dependence of the critical temperature of wires or thin films on the 
(parallel) magnetic field for several values of magnetic impurity concentration. Note the quali- 
tative features first obtained by Kharitonov and Feigel'man [7] : starting at low concentrations 
of magnetic impurities, the critical temperature decreases monotonically with the applied mag- 
netic field. For larger concentrations, a marked non-monotonicity develops, and for yet larger 
concentrations, a regime is found in which the magnetic field first induces superconductivity 
but ultimately destroys it. The physical picture behind this is the competition mentioned in 
the Introduction: first, by polarizing the magnetic impurities the magnetic field suppresses 
their pair-breaking effect. At yet larger fields, this enhancing tendency saturates, and is then 
overwhelmed by the pair-breaking tendency of the orbital coupling to the magnetic field. 

Results for the critical current density. - To obtain the critical current density j Cl we 
first determine the current density (average over the sample thickness) from the solution of 
the Usadel equation via 

j(u) = 2eN TrDTY,Jle(F? 2 (Lu) [u - *{(A))] ) , (23) 

( 1 )We note that, simplifications associated with the strong spin-orbit scattering assumption and the power 
series expansion in F are only necessary to ease the numerical calculations. Our conclusions are not sensitive 
to these simplifications in the parameter regimes considered in Figs. HI and 121 
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Fig. 1 - Critical temperature vs. (parallel) magnetic field for a range of exchange scattering strengths 
characterized by the dimensionless parameter a = Ti/{kBTcoTn)- The strength for potential scattering 
is characterized by parameter ^/(fcsTcoTtr) = 10000.0, and that for the spin-orbit scattering is by 
h/(kBTcoT so ) = 1000.0; the sample thickness is d = 90.0h/pF, where pf is the Fermi momentum; 
the impurity gyromagnetic ratio is chosen to be g s — 2.0; and the typical scale of the exchange energy 
U2 in Eq. &3bl is taken to be Ef/7.5, where Ef is the Fermi energy. 

Fig. 2 - Critical current vs. (parallel) magnetic field at several values of temperature, with the 
strength of the exchange scattering set to be a = 0.5 (corresponding to the solid line in Fig. 0, and 
all other parameters being the same as used in Fig. Q 



and then maximize j(u) with respect to u. In the previous section, we have seen that, over 
a certain range of magnetic impurity concentrations, Tc displays an upturn with field at 
small fields, but eventually decreases. Not surprisingly, our calculations show that such non- 
monotonic behavior is also reflected in the critical current. 

Perhaps more interestingly, however, we have also found that for small concentrations of 
magnetic impurities, although the critical temperature displays no non-monotonicity with the 
field, the critical current does exhibit non-monotonicity, at least for lower temperatures. This 
phenomenon, which is exemplified in Fig. sets magnetic impurities apart from other de- 
pairing mechanisms. The reason why the critical current shows non-monotonicity more readily 
than the critical temperature does is that the former can be measured at lower temperatures, 
at which the impurities are more strongly polarized by the field. 

Conclusion and outlook. We address the issue of superconductivity, allowing for the 
simultaneous effects of magnetic fields and magnetic impurity scattering, as well as spin-orbit 
impurity scattering. In particular, we investigate the outcome of the two competing roles 
that the magnetic field plays: first as a quencher of magnetic impurity pair-breaking, and 
second as pair-breaker in its own right. Thus, although sufficiently strong magnetic fields 
inevitably destroy superconductivity, the interplay between its two effects can, at lower field- 
strengths, lead to the enhancement of superconductivity, as first predicted by Kharitonov 
and Feigel'man via an analysis of the superconducting transition temperature. In the present 
Letter, we adopt the Eilenberger-Usadel semiclassical approach, and are thus able to recover 
the results of Kharitonov and Feigel'man, which concern the temperature at which the normal 
state becomes unstable with respect to the formation of superconductivity; but we are also 
able to address the properties of the superconducting state itself. In particular, our approach 
allows us to compute the critical current and specifically, its dependence on magnetic field 
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and temperature. 

We have found that any non-monotonicity in the field-dependence of the critical temper- 
ature is always accompanied by the non-monotonicity of the field-dependence of the critical 
current. However, we have also found that for a wide range of physically reasonable values 
of the parameters the critical current exhibits non-monotonic behavior with field at lower 
temperatures, even though there is no such behavior in the critical temperature. 

Especially for small samples, for which thermal fluctuations can smear the transition to 
the superconducting state over a rather broad range of temperatures, the critical current is 
expected to provide a more robust signature of the enhancement of superconductivity, as 
it can be measured at arbitrarily low temperatures. In addition, the critical currents can be 
measured over a range of temperatures, and can thus provide rather stringent tests of any the- 
oretical models. Recent experiments measuring the critical temperatures and critical currents 
of superconducting MoGe and Nb nanowires show behavior consistent with the predictions of 
the present Letter, inasmuch as they display monotonically varying critical temperatures but 
non-monotonically varying critical currents [8] . 
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